function [Nshape,dNdzeta]=GetLineNshape(NdNum,zeta)
Nshape=zeros(NdNum,1);
dNdzeta=zeros(1,NdNum);

Nshape(1)=0.5*(1-zeta);
Nshape(2)=0.5*(1+zeta);

dNdzeta(1)=-0.5;
dNdzeta(2)=0.5;
